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ABSTRACT 

We study the gas content of halos in the early universe using high resolution hydro- 
dynamical simulations. We extract from the simulations and also predict based on 
linear theory the halo mass for which the enclosed baryon fraction equals half of the 
mean cosmic fraction. We find a rough agreement between the simulations and the 
predictions, which suggests that during the high-redshift era before stellar heating, 
the minimum mass needed for a minihalo to keep most of its baryons throughout its 
formation was ~ 3 x 10 4 M . We also carry out a detailed resolution analysis and 
show that in order to determine a halo's gas fraction even to 20% accuracy the halo 
must be resolved into at least 500 dark matter particles. 
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1 INTRODUCTION 

The formation of galaxies is one of the most important re- 
search areas in cosmology. Within the simplified hierarchical 
scenario of a universe governed by a cosmological constant 
and cold dark matter, the densi ty profiles of forming ha- 
los h ave been well characterized (|Navarro. Prenk. fc White! 
1 199 it ). However, the complex processes of gas dynamics, such 
as cooling and heating mechanisms, that are responsible for 
the formation of luminous objects still pose many theoretical 
difficulties. 

Numerical calculations show that the first generation of 
galaxies formed at ve ry high redshifts inside collapsi ng halos 
(starting at z ~ 65; iNaoz, Noter. fc Barkanal l2006t ) . corre- 
sponding to high peaks of the primordial density field. In- 
deed, The Wilkinson Microwave Anisotropy Probe ( WMAP) 
measured a Thomson scatterin g optical depth of r e = 0.09 ± 
0.019 from their 5-year data (|Dunklev et 81.1120081 ). When 
combined with simple analytic perscriptio ns of the growth 
of th e global ionized fraction (e.g., fig. 22 iBarkana fc Loebl 
l200ll : IMesinger. Johnson, fc Haimanll2006l . fig. 9), this mea- 
surement suggests that reionization began at very high red- 
shifts, z > 15. This means that a high enough abundance 
of luminous objects must have existed at that time, since 
these first luminous objects are expe cted to have heated and 
reionized their surroundings (e.g.. IBarkana fc Loebl I200TI : 
IWvithe fc Loebl 120031 : lHaiman fe Holderl 120031 ; ICenN2003l ). 
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The formation of a luminous object inside a halo requires, 
of course, gas to be inside the halo. Even in halos that are 
too small for cooling via atomic hydrogen, i.e., minihalos, 
the gas content can have substantial, and observable, as- 
trophysical effects. In addition to the posibility of hosting 
astro physical sources, minihalos may p r oduce a 21-cm signa- 
ture (IKuhlen. Madau. fc Montgomery! (l20q6l);IShapfro et all 
(12003) : INaoz fc Barkanal l|2008l ) but see iFurlanetto fc Oh 
(2006)), and they can block ionizing radiation and pro- 
duce an o verall delay in the global progress of re i oniza- 
tion (e.g.. IBarkana fc Loebl 120021 : llliev et al.l 12001 120051 : 
iMcQuinn et alj|2007l) 7 Thus, the evolution of the halo gas 
fraction at various epochs of the universe is of prime impor- 
tance, particularly in the early universe. 

The estimation of the gas fraction in simulations 
and semi-analytical models has been extensi vely inves- 
tigat e d and used for various purpos e s (e.g.. | Efstathiou 
19921 : IShapiro. Giroux. fc Babul Il994l: iThoul fc Weinberg 



19961 : iQuinn. Katz. fc Efstathioul 1 19961 : iHui fc Gnedin 
19971: iBromm. Coppi. fc Larson! 1 19991: iGnedinl 12000 
20001: 



Kitayama & Ikcuchi 



Abel, Bryan, & Normar 



2002 : IBromm. Coppi. fcLarsonl 120021: iHellv et all |20Q. 



Bromm fc Loebl 12004 iDiikstra et all |2004 lO'Shea et al 



2005al: IMesinger. Bryan fc Haimanl 120061 ; I Naoz fc Barkana 



20071 : IMesinger. Bryan, fc Haimanl I I2OO8I : iTrenti fc Stiavelli 
20091 ). When investigating this issue with simulations, a 
large volume is needed for an adequate statistical sample 
of halos, but on the other hand it i s critical to maintain 
the pr oper resolution. The results of ISpringel fc Hernquistl 
(2003) showed that in order to determine the mass and 
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merger history of each halo even crudely in simulations, 
each halo must be resolved into 500 particles. In the high 
redshift regime, where halos are small and rare, these reso- 
lution requirements are not easy to achieve. Nonetheless, it 
is important to do so in order to understand the formation 
of the first objects. 

Consider the various scales involved in the formation of 
non-linear objects containing DM and gas. On large scales 
(small wavenumbers) gravity dominates halo formation and 
gas pressure can be neglected. On small scales, on the other 
hand, the pressure dominates gravity and prevents baryon 
density fluctuations from growing together with the dark 
matter fluctuations. The relative force balan ce at a given 
time can be characterized by the I Jearisl l)l928h scale, which 
is the minimum scale on which a small gas perturbation 
will grow due to gravity overcoming the pressure gradient. 
As long as the Compton scattering of the cosmic microwave 
background (CMB) on the residual free electrons after cos- 
mic recombination kept the gas temperature coupled to that 
of the CMB, the Jeans mass was constant in time. However, 
at z ~ 200 the gas temperature decoupled from the CMB 
temperature and the Jeans scale began to decrease with time 
as the gas cooled adiabatically. Any overdensity on a scale 
more massive than the Jeans mass at a given time can be- 
gin to collapse, due to a lack of sufficient pressure. However, 
the Jeans mass is related only to the evolution of perturba- 
tions at a given time. When the Jeans mass itself varies with 
time, the overall suppression of the growth of perturbations 
dep ends on a time-averag ed Jeans mass. 

iGnedin fc Hul CL998) 

defined a "filtering mass" that de- 
scribes the highest mass scale on which the baryonic pressure 
still mana ges to suppress the linear baryonic fluctuations sig- 
nificantly. iGnedinl (2000) suggested, based on a simulation, 
that the filtering mass also describes the largest halo mass 
whose gas content is significantly suppressed compared to 
the cosmic baryon fraction. The latter mass scale, in gen- 
eral termed the "characteristic mass" , is defined as the halo 
mass for which the enclosed baryon fraction equals half the 
mean cosmic fraction. Thus, the characteristic mass distin- 
guishes between gas-rich and gas-poor halos. Many semi- 
analytical models of dwarfs galaxies often use the charac- 
teristic rnass_s£a2e_jn_OTdej^_to^^ 

halos (iBullock, Kraytsov. fc Weinberg! l2000l; I Benson et al.l 



Table 1. Main characteristics of the simulations. 



2002a. b; Somervillel I2002T) Recent ly, iHoeft et all <|2006h and 
Okamoto. Gao. fc Theunsl (|2008T ) showed th at the charac- 
teristi c mass scale does not agree with the IGnedin fc Hul 
(1998) filtering mass in the low-redshift, post-reionization 
regime. 

In this paper we explore the very high-redshift regime 
usi ng three-dimensional hydrod y namic al simulations based 
on iMesinger. Bryan fc Haimanl ((2006) . They investigated 
the effects of a photo-ionizing ultraviolet (UV) flux on the 
collapse and cooling of small halos in relic HII regions at 
high redshift, by varying the strength and duration of a tran- 
sient UV and persistent Lyman-Werner (LW) background. 
We consider two different scenarios presented there (see Sec- 
tion [5] for a description of the simulations). We perform a 
resolution study and place a lower limit on the number of 
particles needed in a simulated halo in order to accurately 
determine the gas fraction in halos (Section[5J • We also com- 
pare the numerical evolutio n of the filtering mass a s given 
by the linear calculation in iNaoz fc Barkanal (|2007l ) to the 



Simulation 


Box size 


Static refinement 


Inner refined 




name 


[Mpc/h] 


levels 


region [Mpc/h] 


[M Q ] 


VLres 


1 


2 


0.25 


747 


Lrcs 


1 


3 


0.25 


93 


Hres 


0.35 


3 


0.0875 


4.0 


VHres 


0.25 


3 


0.0625 


1.46 



simulation results (Section |4|. Finally, we summarize and 
discuss our results in Section [S] 



2 THE SIMULATIONS 

Our simulations assume the following cosmological param- 
eters: (O a , (lu, fib, n, <r 8 , H ) = (0.7, 0.3, 0.047, 1, 
0.92, 70 km s" 1 Mpc" 1 ). We use the Eulerian adaptive 
mesh refinemen t (AMR) cod e Enzo, which is describe d in 
greater detail in Brvanl (Il999l ). Norman fc Br yan (1999]) and 
IMesinger. Bryan fc Haimanl 1 20061 ) . Our simulation volumes 
are initialized a t Zj n jt = 99, wi t h den sity perturbations 
drawn from the lEisenstein fc Hul (|l999l ) power spectrum. 
The initial density fluctuations are assumed equal for the 
baryons and dark matter, and the initial gas temperature is 
uniformly 110.7 K. The parameters of our runs (identified 
in order of increasing mass resolution as Very Low, Low, 
High, and Very High resolution) are shown in Table [T] In 
each case, we first run a low resolution dark matter only 
run, in which the simulation volume Zj? oot is resolved on a 



128 root grid. We include an additional n re f static levels of 
refinement (listed in Table [T] for each run) inside this cen- 
tral cube. Furthermore, grid cells inside the central region 
are allowed to dynamically refine so that the Jeans length 
is resolved by at least 4 grid zones and no grid cell contains 
more than 4 times the initial gas mass element. Each addi- 
tional grid level refines the mesh length of the parent grid 
cell by a factor of 2. We allow for a maximum of 11 levels 
of refinement inside the refined central region, granting us 
a spatial resolution of Z root /(128 x 2 11 ). The dark matter 
particle mass for each run is shown in Table [T] We also in- 
clude the non-equilibrium reaction network of nine chemical 



species (H, H^ 



, He, H e" 1 



H e ++ , e~, H 2 , H+, H~) using the 



algorithm of|A nninos et all (Il997|). and initialized with post - 
recombination abundances from lAnninos fc Normanl (fl996). 
Our analysis below is based on the central refined region; the 
low resolution dark matter outside the refined region serves 
to provide the necessary tidal forces to our refined region. 

Each simulation includes two ionization scenarios: (1) 
no UV background radiation (hereafter "NoUV"), and (2) 
flash ionization at z — 25 ("Flash") which instantaneously 
sets the gas temperature to T=10,000 K and the hydrogen 
neutral fraction to xhi = 10~ 3 throughout the simulation 

volume, but involves no heating thereafter. 

We use the HOP algorithm j|Eisenstein fc Hut|[l998l) on 
the DM particles to identify DM halos. Since we are inter- 
ested in the gas fraction, it is essential to sum up both the 
DM mass and the gas mass over the same volume. We obtain 
the total DM and baryonic mass of each halo by integrating 
the densities over a sphere whose radius is the halo's virial 
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radiufl Additionally, we discount halos which have been 
substantially contaminated by the large (low-resolution) DM 
particles outside of our refined region. Specifically, we re- 
move from our analysis halos with an average DM particle 
mass greater than ~ 130% of the refined region's DM mass 
resolution. 

The simulations do not attempt to model global 
reionization. Instead we focus on a short-lived UV back- 
ground that can occur from a localized star forma- 
tion episode. Specifically, the Flash scenario enabled 
iMesinger. Bryan fc Haiman (2006) to compare their results 
to lO'Shea et all l|2005ah ~ and to check the importance of 
including additional dynamical effects of prolonged photo- 
heating. For us, this is simply a good test case of a time- 
variable Jeans mass at high redshift. 



3 THE RESOLUTION DEPENDENCE OF THE 
GAS FRACTION IN HALOS 

Our simulation runs at four different resolutions (Table [TJ 
allow us to carefully test the dependence of the gas fraction 
estimates on the resolution. We find that the halo gas frac- 
tion indeed varies strongly with the resolution. This can be 
seen in Figs. [T] and [2] which show the gas fraction in halos 
as a function of the halo mass. The gas fraction is shown 
for each of our runs, averaged over bins of halo mass. These 
figures show that the lower resolution runs underestimate 
the gas fraction at a given halo mass, often substantially. 
However, our two highest-resolution runs, Hres and VHres, 
agree in their gas fractions over essentially their entire halo 
mass range, and above that range (i.e., at M > 1O 6 M0) the 
two lower resolution runs agree. This demonstrates numer- 
ical convergence and suggests that if we take our highest- 
resolution run at each halo mass, then we obtain the correct 
value of the mean halo gas fraction. 

Getting a correct gas fraction in simulations is im- 
portant in calculating, for example, the galaxy abun- 
dance. Underestimation of the gas fraction in halos 
can give incorrect results. The V Lres run was used by 
IMesinger. Bryan fc Haimanl (|2006l ) and at the time was 
considered a very high resolution run in the high-redshift 
regime. We find, however, that our high resolution runs (i.e., 
Hres and VHres) estimate the gas fraction in the NoUV case 
to be higher than for the VLres run by a factor of ~ 1.6 for a 
halo m ass of 10 5 Mq. Note, that IMesinger. Bryan fc Haimanl 
(2006) focus on whether the gas has cooled relative to other 
runs at the same resolution, not the exact amount of gas 
inside each halo. Furthermore, that study focused on halos 
whose gas was capable of cooling via the molecular hydrogen 
channel, corresponding to masses greater than ~5x 1O 5 M0 
at these redshifts, where the gas fraction is beginning to 
converge to the higher-resolution values. 

We can further quantify the resolution dependence of 



1 Wc define the virial radius as the radius of the spherical volume 
within which the mean density is A c times t he critical density. 
A e is obtained through the fitting formula in iBrvan fe Normaril 
J 1998ft and is ~ 178 in an Einstein-DeSitter universe. We find that 
the halo masses thus obtained generally agree within a factor 
of two with masses obtained directly from the HOP algorithm 
assuming a mean baryon fraction. 
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Figure 1. NoUV case: The halo gas fraction f g as a function of 
the halo mass in the simulations. We consider the following 
simulations runs: VLres (squares), Lres (triangles), Hres (crosses), 
and VHres (circles). Note that we did not use the Lres data at 
z = 20.5 because of a numerical error. 



the halo gas fraction. At each halo mass, we adopt the 
gas fraction as given by our highest-resolution simulation 
available at that mass scale (/ g (mghsst-res)) to be the cor- 
rect value, and compare to this the result from all lower- 
resolution simulations that have halos of that same mass 
at the same redshift. We plot in Figure [3] these ratios, 
Sal ^/g(Highcst-rcs), where each simulation at each output red- 
shift is compared to the highest-resolution run at the same 
halo mass and redshift. In order to test the idea that the 
resolution effect depends primarily on the halo mass reso- 
lution, these ratios are shown as a function not of the halo 
mass but of the number of particles in each halo (Nh) in 
the lower-resolution run. Note that for this we divided our 
various runs into the same mass bins. Now, if the halo gas 
fraction did not depend on resolution then we would simply 
get unity. Instead, the figure shows that the lower resolution 
runs underestimate the gas fractions in halos. In Figure[3]we 
show together all the various redshift outputs from all the 
different runs; all the redshift are consistent with a single 
relation, i.e., they agree quantitatively. Despite the scatter, 
there is a rather uniform trend among all the points (exclud- 
ing the open circles in the Flash case - see below) , confirming 
the idea that the resolution dependence is mainly an effect 
of halo mass resolution. 

The gas fractions in some halos are less accurately deter- 
mined in the Flash case, for a simple reason. In this scenario 
the global heating in the simulation evaporates the gas, i.e., 
the heating raises the characteristic mass (see also sectionU). 
This is easily seen in Fig. [2] where we observe a declining 
gas fraction with time in the low-mass range (consistent with 
iHaiman. Rees fc Loeblll997| ). Thus, the gas fraction in each 
halo is naturally more sensitive to the halo's surroundings 
and to numerical errors when the gas fraction is very low. 
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Figure 2. Flash case: The halo gas fraction f g as a function of the 
halo mass in the simulations. Same conventions as in Fig.fi] 
Note that we did not use the Hres data at z = 20.5 because of a 
numerical error. 
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Figure 3. The ratio between the halo gas fraction in one simula- 
tion and the gas fraction in the highest-resolution run in the same 
mass bin, as a function of the number of particles per halo (iVfc) 
in the lower-resolution run. We show these ratios from all the out- 
put redshifts. We consider the NoUV (bottom panel) and Flash 
(upper panel) cases, with symbols corresponding to the lower- 
resolution run using the same conventions as in Fig. [T] For the 
Flash case we only include halos for which f g > f g>rn i n = 0.03 in 
the higher-resolution run (see text), and we separately show the 
results from the excluded halos (open circles). 



Indeed, we find large scatter below f g ~ 0.03, which cor- 
responds to low-mass halos (below M ~ 5 x 10 4 Mq) at 
redshifts well after the UV flash. We thus consider only the 
halos for which the gas fraction is larger than the fiducial 
value f g ,m,in = 0.03|j The upper panel of Fig. [3] confirms 
the need for this separation, showing that the Flash halos 
with f g > fg.min follow a similar trend as the halos in the 
NoUV case, while the excluded halos with f g < f g , m in are 
inconsistent and show a much larger scatter. 

To derive a useful result from Fig. [3] we first condense it 
efficiently by putting all the points together into a single set 
of bins in Nh- Specifically, we show in Fig. [4] one set of bins 
for the NoUV case, and one for the Flash case (including 
only the halos with f g > f g , m in)- This figure shows a clear 
trend, consistent between the two cases, of an artificially 
declining gas fraction in poorly resolved halos. If we desire 
a reasonably accurate determination of the mean halo gas 
fraction, with a systematic error of no more than 20% (i.e., 
/g// 9 (Hi g hest-rcs) ^ 0.8), then at least 500 particles per halo 
are required. For a 10% error, ~ 2000 particles are required 
per halo. On the other hand, halos resolved into only 100 
particles underestimate the enclosed gas fraction by more 
than a factor of two. We adopt 500 as a minimum number, 
and henceforth consider only halos for which Nh ^ 500. 

We note that, averaged over the output redshifts, the 
number of halos that we extract from each simulation is 104, 
1532, 2479, and 2304, for the NoUV case and 102, 1641, 
2601, and 324 for the Flash case, for the VLres, Lres, Hres, 
and VHres simulation, respectively. We also note that at 
redshift 20.5 we did not use the data for the simulation Lres 
for the NoUV case and Hres for the Flash case, due to nu- 
merical problems in these particular outputs. It is important 
to point out that we only performed this resolution test with 
Enzo. Although gas properties in dense regions are found to 
agree well in AMR and SPH cosmological simulations, dif- 
ferent gravity solvers can lead to non-negligi ble differences 
in the N-body mass function at low masses (O'S hea et al.l 
l2005al ). Specifically, it seems that the gravitational soften- 
ing in Enzo causes small halos to form somewhat late, but 
detailed resolution tests show that the halo gas fractions at 
a given time are unaffected by sof tening at least down to 
100-particle halos (see Figure 14 in lO'Shea et al.ll2005bl ). 



4 THE CHARACTERISTIC MASS AT HIGH 
REDSHIFT 

4.1 Definition and relation to linear model 

In linear theory t he filtering mass, first defined by 
iGnedin fc Hul { 1998), describes the highest mass scale on 
which the baryon density fluctuations are suppressed sig- 
nificantly compa r ed to the dark matter fluctuations. In 
iNaoz fc Barkanl J2007T I we included the fact that the 
baryons have smoother initial c onditions than the dark mat- 
ter (see INaoz fc Barkanall2005l ) and found a lower value of 
the filtering mas s (by a factor of 3—10, de pending on the red- 
shift). Following INaoz fc Barkanal l|2007l ). the filtering scale 



2 We make this separation only in this resolution analysis, while 
in section [4] we calculate the scatter in each mass bin. 
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Figure 4. The ratio between the halo gas fraction in one simula- 
tion and the gas fraction in the highest-resolution run in the same 
mass bin, as a function of the number of particles per halo (N^) 
in the lower-resolution run. We consider the NoUV case (squares) 
and the Flash case (circles), where in each case all the points from 
the various runs and output redshifts shown in Fig. [3] have been 
condensed into a single set of bins. In the Flash case we include 
only halos for which f g > f g>m i n = 0.03 in the higher-resolution 
run. 



where /b,o is the gas fraction in the high-mass limit. In this 
function, a higher a causes a sharper transition between 
the high-mass (constant f„) limit and th e low-mass limit 
(assumed to be f g oc M 3 ). Gnedinl (|2000l ) found a good fit 
for a = 1, with a characteristic mass that in fact equaled the 
filtering mass by his definitio n. By our definition in eq. @, 



Dy nis 

the claim from IGnedinl (120001) is t hat M c = 8 x M F 



Recently, 



lOkamoto. Gao. fc Theuns 



Hocf t et al 



(120081) 



(2006) 
used higher 



and 
resolu- 



Gnedin (2000) and showed that 



tion simulations than in 
this claim is incorrect in the low-redshift regime. They 
compa red the charact e ristic mass found in their simulations 
to the iGnedin "fc Hul { 19981 ) filtering mass and found that 
the two values diverge after reionization. Specifically, they 
found that a = 2 and M c (z = 0) ~ 6.5 x 10 9 M©//i, which 
is much lower than the filtering mass (and thus even lower 
compared to 8 x Mf). We caution that unlike our simple 
test cases, the heating in reionization simulations is complex 
and inhomogeneous, and thus the filtering mass cannot be 
directly and precisely defined and computed. Also note that 
the precise quantitative results during and after reionization 
depend on the thermal history of the gas which observa- 
tionally is not well constrained, though results at z = are 
more robust, as the gas in overdensities around halos has 
had more time to "forget" the details of its heating history. 
Nevertheless, by simulati ng a large set of thermal his tories 
in a cosmological setting, iMesinger fc Diikstral (|2008l ) were 
able to draw a general conclusion that the characteristic 
mass towards the end of reionization is likely in all cases to 
be close to the atomic-cooling threshold of ~ 1O 8 M0. 



(specifically, the filtering wavenumber kp) is defined by ex- 
panding the ratio of baryonic to total density fluctuations 
to first order in fc 2 : 



1 — TJT + r LSS 



(1) 



where k is the wavenumber, and 5b and <5tot are the bary- 
onic and the total (i.e., including both baryons and dark 
matter) density fluctuations, respectively. The parameter 
tlss (a negative quantity) describes the relative difference 
between 5b and Stot o n large scales (for more details see 
iNaoz fc Barkanall2007l ). The filtering mass is defined from 
kF simply as: 

4tt _ (1 2tt\ 3 

M ' = T po \2k;) ' (2) 

where po is the mean matter de nsity today. T his relation 
is one eighth of the definition in IGnedinl ([20001 ) (based on 
a non-sta ndard definition of the Jeans mass used there). 
Following INaoz fc Barkanal i|2007l) we calculate the filtering 
mass for the cosmological parameters assumed in this paper. 

There is no apriori reason to think that the filtering 
mass can also accurately describe p roperties of hig hly non- 
linear, virialized objects. For halos, IGnedinl (|200ol ) defined 
a characteristic mass M c for which a halo contains half the 
mean cosmic baryon fraction /t,. In his simulation he found 
the mean gas fraction in halos of a given total mass M, and 
fitted the simulation results to the following formula: 



/b,0 



1+2 



a/3 



-3/a 



(3) 



4.2 Simulation results 

In order to determine the characteristic mass from our simu- 
lations, we put together the gas fraction measurements from 
all the simulation runs at each redshift, but including only 
the well-resolved halos, i.e., those with Nh 500, as deter- 
mined in section [3] We fit the simulation results to equa- 
tion ([3]) with two free parameters, M c and a, taking /b,o to 
be the average gas fraction in the highest few mass binijf]. 
We used a minimum-x 2 method to estimate the best fit to 
equation (J3j and to find the errors. To account for the nu- 
merical scatter we fitted this equation to the binned data, 
adopting the standard deviation within the mass bin as the 
uncertainty in each binned value. The best-fit parameters 
for both the NoUV and the Flash scenarios along with their 
l-(j confidence limits are listed in Table [2] (and also shown 
in Fig. [7] which is discussed below) . We also compare the 
binned data to the corresponding best fits in the form of 
eq. pt in Figure [5] 



iHoeft et al] l|2006h and lOkamoto. Gao. fc Theunsl 
(2008) found that the fits to their simulations were con- 



3 We note that we also tried to carry out the analysis with fb o 
assumed to be the cosmic mean value. This gave a substantially 
worse fit, with the parameters changed significantly (~ 40% and 
~ 10% changes in M c and a, respectively). We also tried the 
approach of taking fb q to be a free parameter (in the NoUV 
case). This made less of a difference compared to taking fb o 
from the highest few mass bins. We found a lower M c by < 10%, 
a lower a by < 6%, and a fitted value for fbfi higher by < 7% 
than the simple estimate we used in the text. 
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Table 2. The best-fit parameters from equation (|3jl- 



10 



Redshift 


M c 


a 


reduced x 2 degrees 




[10 3 M Q ] 




of freedom 


NoUV 



20.5 

21.5 

22 

22.5 

23 

24 

25 



4.52 
3.28 
5.65 
4.43 
4.83 
3.13 
4.12 



1.45 

1.88 
+0.79 
0.99 
+0.17 
0.47 
+0.64 
0.93 
+0.16 
0.16 
+0.49 
0.47 
+0.69 
1.42 



u - 4 '-0.07 
,+0.03 
-0.03 

70+ 01 

u -' u -0.03 

"■ o,i 0.05 
,+0.02 
-0.02 

75+ ' 21 
U -'°-0.11 
„+0.06 
-0.1 



0.4U 



0.63Z 



0.64"' 



1.15 
0.81 

1.3 
0.54 

1.4 
0.45 
0.21 



32 
52 
49 
46 
35 
33 
21 



Flash 



20.5 


141.3+-° 


0.49 


21.5 


68.7+.°^ 


0.64 


22 


51.47+j;^ 


0.48 


22.5 


55.99 + °23 


0.63 


23 


22.83+J^ 


0.62 


24 


16 2+ 34 


0.72 


25 


4 12+069 
^■ lz -1.42 


0.64 



+0.02 
-0.02 
+0.01 
-0.01 
+0.02 
-0.02 
+0.01 
-0.01 
+0.03 
-0.03 
+0.07 
-0.07 
+0.06 
-0.1 



0.9 

1.6 
3.8 
1.58 
0.85 
0.96 
0.17 



31 
48 
54 
53 
36 
31 
26 
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Figure 5. The gas fraction in all the simulations vs. the halo 
mass at various redshifts. The simulation data have been binned 
by halo mass. We use only data points from halos with ^ 500. 
We consider the Flash (crosses) and NoUV (filled squares) cases. 
Also shown are the best fits in each case in the form of eq. llBl . We 
note that the bins show here are equally spaced for representations 
reasons, for the actual degrees of freedom see table [2] 



sistent with a = 2. However, we find that our fits yield 
Oi ~ 0.4 — 0.7 (see Tableland also Fig. [7] middle panel). In 
Figure |S] we show the dependence of the reduced \ 2 on the 
fit parameters M c and a at z = 20.5. This figure shows that 
a low a gives the best fit to equation (3), and suggests that 



-o 

9) 
o 

3 



0.1 



0.5 




5x10 s 



2xl0 4 



Figure 6. The dependence of the reduced \ 2 on the parameters 
of the fit at z = 20.5. There are 33 degrees of freedom in this case. 
Bottom panel: dependence on M c obtained while fixing a = 2, 1, 
or 0.47, from top to bottom respectively, where the lowest value 
is the best-fit one (see table Upper panel: dependence on a, 
fixing the best-fit value M c = 4.52 X 10 3 M . 



the characteristic mass scale found assuming a — 2 would 
be an underestimate by about 15% in this case. 

Figure [7] shows the best fitted parameters at various 
redshifts for M c and a, and our value for fb,o, for both the 
NoUV and Flash scenarios. The 1-a (68%) confidence re- 
gions are listed in table [5] since these statistical errors are 
small, we do not show them in the figure. Indeed, particu- 
larly in the Flash case, the fits do not drop as quickly as the 
data points do towards low halo masses, and so tend to sys- 
tematically underestimate the characteristic mass. Thus, to 
obtain a more realistic estimate of the systematic error re- 
sulting from the form of the fit, we also show M c as derived 
directly from the binned data without a fit; in this case we 
simply found the maximum halo mass for which f g = /t,,o/2 
(interpolating inbetween the binned points). We find that 
the systematic errors for the NoUV case (typical difference 
of a factor of 1.5-2 between the fit and the no-fit values) are 
substantially smaller than in the Flash case (typical differ- 
ence of a factor of 3), and they are much larger than the 
statistical errors in both cases. 

Figure [7] also shows the analytical calculation from lin- 
ear theory of th e filter ing mass, as described in detail in 
iNaoz fc Barkanal (|2007l ); we make this calculation for both 
the NoUV and Flash scenarios, assuming the same initial 
conditions as in the simulations. As mentioned above, the 
simulation assumes equal baryon and dark matter fluctu- 
ations at its Zi n it = 99 (as is commonly assumed in the 
literature), while the correct baryonic initial conditions are 
smoother (see below). Note that due to the simplicity of the 
Flash scenario that we have implemented in the simulation, 
we can easily incorporate it precisely within the analytical 
calculation. We also directly tested the effect of radiative 
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Figure 7. The parameters of the best fits in the form of equa- 
tion \3\ different panels show M c , a, and /b o- We consider the 
Flash (crosses) and NoUV (filled squares) cases, where we fit 
equation \3\ to all data points from halos with 500. For 

comparison we also show the M c values derived directly from in- 
terpolating the binned data, without assuming a fitting function 
(dotted curves in the bottom panel; bottom: NoUV, top: Flash 
scenario). In the bo ttom panel we als o show the analytical calcu- 
lation according to iNaoz fc Barkanal (120071 1. assuming the same 
initial conditions as in the simulations (short-dashed curves; bot- 
tom: NoUV, top: Flas h). The full calculation assuming the true 
initial conditions as in lNaoz fc Barkanal J2007T> is also shown for 
the NoUV case (long-dashed curve). We note that the NoUV and 
Flash cases mostly overlap in the top panel, which also shows the 
cosmic mean baryon fraction (horizontal dotted line). 



cooling, which is included in the simulations but not the an- 
alytical model. We carried out an additional Flash run with 
a resolution identical to that of the Hres run, but where 
the radiative cooling was eliminated, leaving only adiabatic 
cooling and Compton heating. We found that the gas frac- 
tions (and thus the fitted characteristic mass) did not change 
significantly, and thus verified that radiative cooling has a 
negligible effect on our results. 

We find that, given our systematic errors, the filtering 
mass from linear theory is consistent with the characteris- 
tic mass fitted from the simulations for both the NoUV and 
the Flash casefl It is important to emphasize that in this 
statement we are referring to our definition in equation ©, 
which is one eighth of the original definition which iGnedinl 
(2000) claimed was a good fit to the characteristic mass. 
In any case, we conclude that at least in a particular red- 
shift range (z — 20 — 25) the filtering mass provides a fairly 
good estimate to the characteristic mass, either before stel- 
lar heating or in its initial stages. Since we have not probed 
a larger range of redshifts, we cannot generalize this conclu- 



4 We note that equation (|3j is successful in giving a reduced x2 
of order unity in all cases except one (see Table [2} ■ 



sion. Also, the large systematic errors (particularly in the 
Flash case) reduce the significance of the above conclusion. 

In Figure [7] we also show the filter ing mass from the 
full calculation of lNaoz fc Barkanal ( 2007 ) assuming the cor - 
rect initial conditions (see figure 1 of lNaoz fc Barkanall2005l ), 
in the NoUV case. The correct initial conditions cannot be 
fully directly incorporated in a simulation without starting 
at much higher redshifts than simulators are used to. In 
particular, these initial conditions include the fact that at 
z — 1200 the baryons are still essentially uniform (on scales 
relevant for galactic halos) due to their just ended strong 
coupling to the CMB photons. Based on the agreement we 
have found between the linear theory and the simulations 
for the case of the simulations' initial conditions, we suggest 
that we can estimate the real characteristic mass in the uni- 
verse based on our analytical filtering mass calculation with 
the true initial conditions. 



There are several differences between our simulations' 
initial conditions and the true ones. The initial conditions in 
the simulations assumed a lower temperature (by ~ 30% at 
Zinit = 99) than the exact calculation with Compton heat- 
ing, resulting in lower gas pressure and thus a lower filtering 
mass than in the full analytical calculation. The assump- 
tion that the baryon perturbations follow the dark matter at 
Zinit = 99 creates tendencies to both raise the filtering mass 
(since the baryon fraction, and thus the gas pressure, is too 
high within perturbations) and lower it (since the filtering 
mass reflectes the integrated effect of pressure, and the in- 
tegral is only begun at Zi n a = 99 instead of at z = 1200). 
Also, we note that the best-known cosmological parameters 
are slightly different from those in our simulation, i.e ., in ou r 
analytical calculation we use those of ISpergel et al. | (|2007t ): 
(f2 A , SIm, fib, n, <7 8 , H ) = (0.701, 0.299, 0.0478, 0.957, 0.82, 
68.7 km s" 1 Mpc" 1 ). However, changing the cosmological 
parameters does not play an important role. This is because 
the definition of the filtering mass is independent of as which 
is the major difference in the cosmological parameters. Thus, 
we conclude that the true minimum mass needed for a halo 
to keep at least half of its baryons, in the era before stel- 
lar heating (i.e., corresponding to the NoUV case), is about 
2.7 x 10 4 M at z ~ 20 - 25. 



Finally, we can also use our simulations to look beyond 
the tight fg(M) relation of equation @, and cons ider the 
distrib ution of gas fractions for a given halo mass. I Gnedinl 
(2000) showed that this distribution in the simulation is well 
approximated as a lognormal distribution. Performing the 
same analysis for the NoUV scenario, we find the same, as 
shown in Figure [8] Since we have a limited number of halos, 
we collect all of our data at each redshift and consider the 
ratio between the measured f g and that predicted by equa- 
tion ([3]) with the best-fit parameters. Thus, we assume that 
this relative distribution does not vary strongly with halo 
mass. Since it is easier to deal with a normal distribution, 
we plot the distribution of ln(/ ff // 9jCa ic) and fit to it a nor- 
mal distribution with mean fi and standard deviation a. We 
show the best-fit parameter values with their 1-a confidence 
ranges in Table [3] 
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ln(f« / f,.cJ 

Figure 8. The distribution of gas fractions with respect to the 
prediction of equation J3t . The histograms show the binned data 
of ln(/g// g , ca ic), where / SiCa lc f° r each halo mass is taken from 
eq. (O assuming the best-fit parameters as given in Table [2] We 
also show in each case the best fit to a normal distribution (solid 
curves). 



Table 3. The parameters of the best-fit normal distributions to 
the histograms shown in Figure [8] 



redshift 




(7 


reduced \ 2 


20.5 


0.0003 ±0.015 


0.24 ±0.01 


0.98 


21.5 


-0.14 ±0.014 


0.22 ±0.01 


0.99 


22 


-0.12 ±0.01 


0.22 ±0.01 


0.97 


22.5 


-0.14 ±0.01 


0.21 ±0.006 


0.99 


23 


-0.12 ±0.01 


0.24 ±0.009 


0.98 


24 


-0.20 ±0.02 


0.24 ±0.01 


0.96 


25 


-0.13 ±0.02 


0.22 ±0.01 


0.93 



5 CONCLUSIONS 

We have used three-dimensional hydrodynamic simulations 
to investigate the resolution requirements needed to deter- 
mine correctly the gas fraction in halos in the early uni- 
verse. We considered both a NoUV case with no stellar heat- 
ing, and a Flash case with instantaneous stellar heating. We 
found that the gas fraction in halos is strongly dependent 
on the mass resolution of the simulation (see Figs. [1] and [2} 
both in the NoUV and Flash cases. Using our multiple runs 
at various resolutions, we demonstrated convergence in the 
estimated gas fractions over a wide range of halo masses; 
thus we concluded that these estimates are likely correct. 
Comparing these converged values to the results from lower- 
resolution simulations, we showed (see Figs. [3] and [4| that 
halos that are poorly resolved (in terms of the number of 
dark matter particles) yield artificially low gas fractions. In 
particular, we concluded that to ensure a gas fraction that 



is unbiased to < 20%, there must be at least 500 particles in 
each halo (and 2000 particles for < 10% bias). We showed 
that such a simple condition is a consistent description of 
the resolution dependence over the full range of redshifts 
and heating stages that we investigated. 

We found from the simulations the characteristic mass 
scale below which a halo does not contain most of its baryons 
and can be considered "gas poor". Specifically, we fitted 
eq. (|3]) to the data from the simulations with two parame- 
ters, the characteristic mass M c and a. We found that they 
are roughly constant with redshift for the NoUV case (see 
Fig. [5] for the fit, also Fig. [7] and table [2] for the best fit 
parameters). We compared between the simulations' char- 
acteristic mass and the analytical filtering mass, using the 
simulations (with their particular initial conditions) as case 
studies. We found that the M c estimations from the sim- 
ulations are consistent with the filtering mass from linear 
theory, according to our definiti on which is one eighth of 
the definition that iGnedinl {2000) claimed was a good fit to 
the characteristic mass. Due to our limited redshift range 
and fairly large systematic errors, we cannot be sure how 
general this consistency may be. We note that this agree- 
ment between the simulations and the linear theory in the 
NoUV and the Flash scenarios occurred in a regime where 
the heating is simple and easily incorporated within linear 
theory. In more complicated situations with inhomogeneous 
heating from astrophysical sources, the filtering mass can- 
not be directly and precisely computed, and this may ex- 
plain ap parent inconsistencies between the theory and simu - 
lations (|Hoeft et al.ll2006l ; lOkamoto. Gao. fc Theunsl 120081 ). 
Note that these authors used the iGnedinP 2000l ) definition 
for the filtering mass, which is eight times our definition. 

The agreement between the linear calculation of the fil- 
tering mass and the simulations suggests that using the cor- 
rect initial conditions we can calculate the true minimum 
mass in the real universe for which a halo keeps most of its 
gas during its formation. We calculated the filtering mass 
using the cor rect initial conditions s tarting from very high 
redshift (see. iNaoz fc Bar kana 2005), finding a characteris- 
tic mass of 2.7 x 10 4 Mq assuming no stellar heating prior 
to z ~ 20. (see Fig. [7]l . We note that the filtering mass for 
the Flash scenario with true initial conditions would have to 
be corrected by approximately the same factor as M c in the 
NoUV case; we did not show this curve in Fig. [7]to avoid a 
busy figure. 

As noted in the introduction, minihalos can have im- 
portant, observable effects, in particular on the early stages 
of cosmic reionization. We have taken an important step to- 
wards understanding the importance of minihalos by estab- 
lishing which ones contain substantial amounts of gas. Also, 
since the minimum mass for molecular hydrogen cooling is 
~ W 5 A1q, somewhat higher than the characteristic mass for 
minihalos, we conclude that the gas fraction within the host 
halos of the first stars could be slightly (though probably 
not greatly) reduced compared to that of more massive ha- 
los; this effect is likely missing or quantitatively inaccurate 
in many simulations of the first stars, due to their inabil- 
ity to start with the correct initial conditions at sufficiently 
high initial redshifts. 
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